Generating Entangled Microwave Radiation Over Two Transmission Lines 
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Using a superconducting circuit, the Josephson mixer, we demonstrate the first ex- 
perimental realization of spatially separated two-mode squeezed states of microwave 
light. Driven by a pump tone, a first Josephson mixer generates, out of quantum vac- 
uum, a pair of entangled fields at different frequencies on separate transmission lines. 
A second mixer, driven by a 7r-phase shifted copy of the first pump tone, recombines 
and disentangles the two fields. The resulting output noise level is measured to be 
lower than for vacuum state at the input of the second mixer, an unambiguous proof of 
entanglement. Moreover, the output noise level provides a direct, quantitative measure 
of entanglement, leading here to the demonstration of 6 Mebit.s -1 (Mega entangled 
bits per second) generated by the first mixer. 

PACS numbers: 



Pairs of entangled electromagnetic fields propagating 
on physically separated channels constitute an essential 
resource in quantum information processing, communi- 
cation and measurements [1. 2J. They can be realized 
by squeezing a vacuum state shared by two spatially 
separated modes. This entanglement is revealed in the 
cross-correlations between well chosen quadratures of the 
two fields which fall below the level of quantum vacuum 
noise. Given the considerable development of microwave 
quantum optics, these Einstein-Podolsky-Rosen (EPR) 
states, or spatially separated two- mode squeezed vacuum 
states, have become highly desirable at such frequencies. 
At optical frequencies, EPR states are usually prepared 
by parametric down-conversion of a pump tone using a 
nonlinear medium [H H]. At microwave frequen- 
cies, only single-mode squeezing and two-mode squeez- 
ing between sidebands of a single transmission line have 
been demonstrated so far, using degenerate Josephson 
parametric amplifiers jSHH]- Recently, a dissipationless, 
nondegenerate, three- wave mixer for microwave signals 
based on Josephson junctions was developed |10ffl2] (see 
Figs. 1,2). Strong correlations between the spontaneously 
emitted radiations from two ports have been observed in 
the parametric down-conversion mode [T3], but the ex- 
periment did not prove directly the presence of entangle- 
ment in the separated output fields. 

Here, we describe an interference experiment demon- 
strating that nondegenerate Josephson mixers can entan- 
gle and disentangle usable EPR states of microwave light 
(Fig. 1). A first mixer, called the "entangler", is driven 
by a pump tone while its two input ports are terminated 
by cold loads ensuring that only vacuum quantum noise 
enters the device. The two entangled output ports feed 
the input ports of a second mixer called the "analyzer". 
The role of the analyzer is to recombine and disentan- 



gle the two microwave fields before sending them to a 
standard microwave amplification and detection chain. 
As the phase difference between both pumps varies, the 
noise at the output of the analyzer exhibits interference 
fringes which pass under the level of amplified vacuum. 
Remarkably, the measurement of the noise at the output 
of the analyzer directly quantifies entanglement between 
its two input fields without resorting to two homodyne 
detection channels and the analysis of their correlations. 

The Josephson mixer |10f[T2] is a superconducting cir- 
cuit parametrically coupling two superconducting res- 
onators (Fig. 2) at distinct frequencies f a and ft, via 
a pump at their sum frequency fp = f a + fb- Each res- 
onator has only one access port, but input and output 
signals are spatially separated by cryogenic microwave 
circulators (Fig. 2 and |14j ) so that the entangler out- 
put can be exclusively sent to the analyzer input. Each 
mixer performs a reversible transform of the wavefunc- 
tion of the field via the unitary two-mode squeeze oper- 
ator S = exp {re itpp a)tf - re~ lipp ab) where re iipp is the 
complex squeezing parameter, and a and b are the field 
operators of the two modes |15j . The input and output 
canonical field operators are related by the scattering re- 
lations 

aout = S^a in S = cosh(r)a in + e %H>p smh(r)b\ n , . 
bLt =S^b\ n S= cosh(r)&L + e-^sinh(r)a m ( ' 

where ipp is the phase of the pump and G — cosh 2 r = 
(Pth + P) 2 /(Pth ~ P) 2 is the power direct gain which 
increases with pump power P below the parametric 
self-oscillation threshold P t h- With the pump on, the 
vacuum state at the input is converted into a two- 
mode squeezed vacuum state \Sq) — S\0) \0) b = 
cosh(r) -1 ^tanh(r)™ \n) \n) b . Note that this entangled 
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Figure 1: Principle of the experiment, a. When pumped with a microwave tone at frequency fp~f a + ft, a Josephson 
mixer (black diamond) transforms incoming quantum vacuum noise (left) on modes a and b into an EPR state (right). The field 
states are represented by their standard deviation contours in the single mode phase space of a (top), b (bottom) and in the 
bipartite phase space (middle) spanned by (Re(a), Re(6)) and (Im(a), Im(fc)) where Re(a) = (a + a" ) /2 and Im(a) = (a - a f ) /2i 
are the in-phase and out-of-phase quadratures of mode a. In each plot, a solid black circle sets the scale of vacuum noise and a 
dashed circle sets the scale of amplified vacuum noise. Quantum entanglement in the output fields is observed in the bipartite 
phase space where cross-correlations go beyond quantum uncertainty by a squeezing factor e~ r . b. In order to demonstrate 
entanglement at the ouput of this mixer named "entangler", a second, identical mixer named "analyzer" is placed in series and 
pumped by the same tone with phase difference A<p. Entanglement at the input of the analyzer is revealed by measuring, at 
the analyzer output (dashed frame), a lower noise level on mode a or b than for amplified vacuum (dashed circle) for a given 
phase Aip. The smallest (largest) output noise level occurs for opposite (equal) squeezing factor at Aip — tt (0), and can ideally 
go as low as the vacuum noise level. 



state can be understood as the superposition of twin pho- 
tons with different frequencies and propagating on spa- 
tially separated transmission lines. Non-local two- mode 
squeezing directly appears in the combinations of output 
fields 
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which, for <pp — 0, implies cross-correlations between 
Re(a) and Re(fr) on one hand and Im(a) and — Im(6) on 
the other hand, beating the Heisenberg limit of vacuum 
quantum noise, as shown in Fig. 1. In optics, these 
correlations have been observed by double balanced ho- 
modyne detection techniques in several systems |16j . The 
present experiment describes the first demonstration at 
microwave frequencies of these quantum correlations be- 
tween signals on spatially separated transmission lines. 
The Josephson mixer here serves two functions. First, 
the entangler produces EPR states of microwave light at 
incommensurate frequencies and on spatially nondegen- 
erate modes with squeezing parameter te- Second, the 
analyzer recombines input fields as shown in Eq. , with 
squeezing parameter ta and relative pump phase Aip, in 
order to reverse the transformation and disentangle the 
field state (Fig. 1). 

The output noise of the entangler can be measured 
on each mode independently by turning off the analyzer 



{ta = 0) . The noise power spectrum measured by a spec- 
trum analyzer is proportional to the symmetrized vari- 
ance of the field operator [T7] 

({aout,E,aL t , E }) coshjrg 
(Aa outjB ) = \\a>out,E)\ = 2 ' 

(3) 

The variance of this "amplified vacuum" is always larger 
than that of the vacuum state, for which (Aa) 2 = 1/2 
(Fig. 1). Discarding the information from the other 
mode, each output field is in a thermal state [18 j . Yet, 
since the combined two-field state \Sq) is a pure state 
with no entropy, it is possible, ideally, to reverse the 
squeezing with a second mixer and re-obtain a vacuum 
state on each port. The analyzer can perform this in- 
version if operated with opposite squeezing parameter 
r_4 = — Te- In practice, unavoidable losses between the 
two mixers prevent the exact recovery of the vacuum. 

Losses are modeled as field splitters coupling uncorre- 
cted cold thermal baths to each mode (Fig. 3) so that 



,.4 = V aa out ^ E +Vaa t h and b in A = V pb ou t,E+V 
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where a t y t and bth describe bosonic modes of thermal 
baths at frequencies f a and fb, and a = 1 — a, /3 = 
1 — p. Besides, microwaves photons propagate for a finite 
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Figure 2: Schematics of the experimental setup. Each Joseph- 
son mixer consists of a ring of Josephson junctions coupling 
two A/2 superconducting resonators addressed via a 180° hy- 
brid coupler or a single ended port. Both mixers are de- 
signed with the same geometry as in Ref. 12J, and their 
resonance frequencies are matched at f a — 5.578 GHz and 
fb = 8.812 GHz using two independent flux bias. The pump 
frequencies are set to fp — 14.390 GHz. Microwave circu- 
lators separate the input and output of the entangler and 
analyzer. Input ports are represented as open circles, and at 
each output port, the double triangle symbolizes the low noise 
amplifying measurement setup with total gain Glna- 
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Figure 3: a. Protocol of the scattering coefficient measure- 
ments by a vector network analyzer connected between the 
a, b input and the a, b output ports. The setup is calibrated 
by turning on and off each Josephson mixer separately. Losses 
are modeled as field splitters of transparency a 2 and /3 2 cou- 
pling a cold load to the signals, b. Color traces: Transmis- 
sion measurements of |t a ^ a | 2 as a function of phase difference 
Aip between both pump signals. The gain of the analyzer is 
set to Ga = cosh 2 ta = 10 (solid gray line) . Each trace 
and color corresponds to a different gain for the entangler 
G E = cosh 2 r E = 0.2, 0.8, 1.8, 3.2, 5, 7.2, 9.8 dB. Dashed lines: 
fits to the data using equation Q and the single fit parameter 
P/ a — 0.945. Together with an independent, in situ noise cal- 
ibration this value leads to a = 0.33±0.05 and /3 = 0.36±0.05 
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amount of time r a and between the two mixers leading 
to a correction of the phase difference entering the scat- 
tering terms Aip' — Aip — 2irf a T a — 2irfbTb- The temporal 
extent of the twin photons exiting the entangler is given 
by the inverse of the bandwidth Af = Afo/ coshrg |13j . 
In the experiment, the travel times r a and r& of order 
2 ns are much smaller than this temporal extent since 
A/o = 28 MHz, so that microwave photons do interfere 
even if their travel durations may slightly differ between 
modes. It is then straightforward to calculate the scat- 
tering coefficients of the full circuit. For instance, the a 
output mode is given by 

Q>out,A = ta—^a^injE ~t~ ^b— >a^^ n e 

+\fa coshrvtdt/j + e* Av \ff3 sinhrvi&L. 

where 

t a ^a = cosh r e cosh r a + e lAip y/5 sinh te sinh ta , 
tb^a = v^sinhrg cosh ta + e lAv> -\//3 cosh te sinh r^. 

(4) 

These scattering coefficients were measured using a 
nonlinear four-port vector network analyzer as a func- 
tion of the phase difference Aip for various values of the 
gains cosh 2 te,a ranging from 1 to 40, a subset of which is 
shown on Fig. 3. The special cases where one or both of 
the converters are not pumped (r — 0) offer the opportu- 
nity to calibrate each converter gain independently. The 
only fit parameter for this whole set of measurements is 



the ratio between transmissions on both arms, found to 
be /3/a = 0.945. 

In Fig. 3, only mean values of the output field ampli- 
tudes are measured. Yet, truly quantum features appear 
in their correlations. Consider the case of a cold load 
setting the vacuum quantum state at the input of the 
entangler, which is reached in our experiment at 45 mK 
since hfb/k > hf a /k = 260 mK [14 . When the entangler 
is turned off (te — 0), the analyzer is fed by vacuum fluc- 
tuations and the output noise reads the amplified vacuum 
level (Aa ou t,A) 2 = cosh(2r y i)/2 as in Eq. In general, 
the output noise a 2 , normalized to that reference level, 
on both output ports can be calculated from Eq. ^ and 
oscillates with phase Aip as 

a 2 (Aip) ee 2(Aa out!7l ) 2 /cosh(2r A ) 

= a (cosh 2te + sinh 2rg tanh2r^ cos Aip) + a. 

(5) 

For simplicity, this expression is only given in the case 
of balanced losses a — j3, but the general case can be 
treated without much difficulty. The maximal and mini- 
mal values of a 2 , corresponding to the extrema of cross- 
correlation between quadratures of the EPR state (el- 
lipses in Fig. la), are obtained respectively for Aip = 
and Aip = tt 

2 n ,cosh(2r E ±2r A ) , . 

Omax.min = I 1 ~ a ) T7^ \ + a - 6 

cosh(2r j 4j 

The existence of a phase Aip for which the output noise 
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Figure 4: a. Color traces: variance of the output mode 
(AaoutA) 2 referred to the case of vacuum input on the an- 
alyzer (divided by cosh(2rA)/2) as a function of phase differ- 
ence Aip, determined by measuring the spectral density of the 
noise at the analyzer a output when only quantum noise en- 
ters the entangler a,b inputs. An absolute calibration allows 
exact conversion between both quantities with an error of at 
most ±2.5% |14| . Each color corresponds to the same gain of 
the entangler Ge as in Fig 3 with a fixed gain on the analyzer 
Ga = cosh 2 ta = 10. The horizontal line at a 2 — 1 represents 
the measured noise for amplified vacuum at the output of the 
analyzer (te = 0). For Aip close to ir, the measured noise 
goes below this level, an evidence of entanglement. Dashed 
lines: predicted variance using Eq. extended to the un- 
balanced loss case using a = 0.37 and /3 = 0.40. b. Dots: 
Noise level measured at Atp = (anti-squeezing) and Aip = ir 
(squeezing) as a function of gain Ge for Ga = 10. The size of 
the dots is larger than the error bar. Solid lines: prediction 
using Eq. Q, extended to unbalanced losses as in panel 4a. 
Colored area: consistent values of the noise using the uncer- 
tainty in the calibration of the losses a and /3 |14j . Dashed 
lines: same prediction but without losses, a = /3 = 0. c. 
solid dots: logarithmic negativity measure of entanglement, 
with errors bars, solid squares: entanglement purity. Lines: 
theoretical predictions using the parameters of panel b. 



a 2 goes below 1, which demonstrate correlations beyond 
quantum uncertainty, is a sufficient evidence of entangle- 
ment [HI [131120]. 

The normalized noise power a 2 is obtained by measur- 
ing the spectral density S a (detailed in supplementary 
information |14| ) 



a 2 (Aip) 



S a (Aip) ~ SoS . 1 



(7) 



cosh(2r^) V hf a G LNA 2 J ' 
where the noise background due to the following ampli 



fiers S s is small enough to be precisely subtracted. The 
spectral densities S a (f a ) and Sb(fb) of both modes at 
the output of the analyzer were measured using a mi- 
crowave spectrum analyzer behind a cryogenic low-noise 
preamplifier on a 0.5 MHz bandwidth. This bandwidth 
was chosen to be smaller than that of the mixers for all 
combinations of gains Ge.a — cosh re, 4 and phase dif- 
ferences Aip. Importantly, it was possible to calibrate 
the total gain of the measurement setup Glna so that 
the normalized noise power a 2 is measured with at most 
±2.5% relative error |14| . This calibration was performed 
by turning on a single mixer at a time and varying the 
temperature of a thermally decoupled input load on mode 
din,E- As a side result, the calibration provides the loss 
a = 0.33 ± 0.05 between mixers on mode a which, to- 
gether with the ratio ft/a — 0.945 from Fig. 3, leaves no 
unknown parameters in the experiment. 

As can be seen on Fig. 4a, the noise does pass below 
the threshold of amplified vacuum noise, hence proving 
the existence of entanglement. Note that measurements 
on mode b (not shown) gave similar results as expected. 
Note also that minimum noise a 2 ^ occurs at |rg| < \ta\ 
and not at exactly opposite squeezing parameters. This 
deviation may be due to the beginning of a saturation 
of the analyzer mixer, corroborated by the slight devi- 
ations of the fits in Fig. 3. For each squeezing param- 
eter te, it is possible to extract the extrema of noise 
c 2 ;„ from the curves of Fig. 4a. These extremal 
noise measurements (Fig. 4b) are well described by Eq. 
^ generalized to unbalanced losses between modes with 
a = 0.37 and j3 — 0.40, consistently with the calibration. 
It is even possible to account for the whole dependence 
of the measured noise on phase difference Aip by gen- 
eralizing Eq. (|5| using the same parameters (Fig. 4b). 
The overall minimum for the measured noise is reached 
at cosh 2 rs « 5 and reads o"^ in = 0.45 ± 0.01 with a 
corresponding maximum er 2 lax — 11.9 ± 0.1. 

It is remarkable that the amount of noise at the out- 
put of a single port of the analyzer directly measures 
the entanglement between the two input fields. In par- 
ticular, the minimum of output noise is linked to the 
logarithmic negativity En = — \og 2 {(T 2 ain ) — 1.15 ± 0.04 
and to the entropy of formation Ef — 0.69 ± 0.03 en- 
tangled bits (ebits, see Q3]) fllA |2TH2"5] : the deeper the 
noise fringes, the larger the entanglement. The purity 
of the entangled state is also related to both extrema 
tr(p 2 ) = (al^cr^)- 1 = 0.186 ± 0.09. These quantities 
of entanglement are within a factor of 2 from the state of 
the art in optics [2T| [26H28] . Given the bandwidth of the 
mixers, the analyzer receives a usable rate of 6 Mebits.s -1 
from the entangler. 

In conclusion, we have demonstrated the production 
of EPR states of microwave radiation. Vacuum noise at 
the input of a first mixer is converted into two entan- 
gled fields. A second mixer is used to recombine and 
disentangle the two fields. Using an absolute calibra- 
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tion, the minimal noise intensity at the output of the 
second mixer, when the phase difference A.ip is varied, 
constitutes a direct measure of the entanglement between 
the twin fields. Our measurements arc limited by the fi- 
nite losses between mixers but still show that a rate of 
6 Mebits.s -1 travel between the entangler and the an- 
alyzer. This first implementation of spatially separated 
two-mode squeezed states in the microwave domain opens 
novel experiments on quantum teleportation or super- 
dense coding in the fields of nanomechanical resonators 
and superconducting qubits. Moreover, inserting a "cir- 
cuit QED" readout cavity in one arm of the vacuum 
quantum noise interferometer described in this paper, 
one would achieve a maximally efficient measurement, 
for a given photon number, of the phase shift associated 
with a change of qubit state. 
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